Chiral /-wave Topological Superfluid in Triangular Optical Lattices 
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We demonstrate that an exotically chiral /-wave topological superfluid can be induced in cold- 
fermionic-atom triangular optical lattices through the laser-field-generated effective non-Abelian 
gauge field, controllable Zeeman fields and s-wave Feshbach resonance. We find that the chiral 
/-wave topological superfluid is characterized by three gap less Majorana edge states located on the 
boundary of the system. More interestingly, these Majorana edge states degenerate into one Majo- 
rana fermion bound to each vortex in the superfluid. Our proposal enlarges topological superfluid 
family and specifies a unique experimentally controllable system to study the Majorana fermion 
physics. 
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Introduction: Topological superconductors/ superflu- 
ids (TSCs/TSFs)[r,'2] have attracted considerable inter- 
est in condensed matter physics because of their potential 
applications on the fault-tolerant topological quantum 
computation (TQC)[3]. One of the remarkable features 
in TSCs/TSFs is the helicity or chirality of the uncon- 
ventional pairings. Unfortunately, there are very limited 
natural materials [3, ^ exhibiting these kinds of uncon- 
ventional pairings. Starting with the pioneer work by Fu 
et al. y, recently, some new classes of hybridized sys- 
tems jTHlOj have been proposed as possible candidates 
for TSCs, where the unconventional pairings are induced 
from proximity effects of conventional 5-wave supercon- 
ductor films. However, the impurities or disorders in the 
materials hosting the electron gas increase the difficulties 
to investigate the topological properties in the hybridized 
systems from experiments. Therefore, it should be not 
only interesting but also necessary to design other sys- 
tems that present TSC/TSF phases. 

The ultracold atom gas associated with optical lattice 
technology provide an ideal platform to realize and in- 
vestigate the topological phases [lll-[l5j due to the con- 
trollability and cleanity. In particular, some chiral p- 
wave TSFs[l2|,|l3 have been proposed based on the laser- 
induced artificial gauge fields [TgI, [I^ in cold atom sys- 
tems. The effect of the artificial gauge fields is equivalent 
to the spin-orbit coupling, a key factor to induce topo- 
logical phase. More recently, some experimental groups 
reported the realization of strong spin-orbit coupling in 
ultracold fermionic atoms gas ^^K and ^Lifis', ^19^. This 
new technique brings the huge hope to realize many ex- 
otic states related with spin-orbit coupling. 

In this paper, we propose that an exotically chiral /- 
wave TSF can be realized through the effective Rashba 
spin-orbit coupling(RSOC) [2Q| . Zeeman field(ZF) and 
5-wave Feshbach resonance in triangular optical lat- 
tices(TOLs). The effective Rashba SOC and ZF are 
produced by the laser-atom interactions through modu- 
lating applied laser beams. The 5- wave Feshbach res- 



onance is utilized to induce the SF states [21] of the 
trapped atoms. We find that there exists a phase 
transition separating the TSF and normal superfluid 
(NSF), which is determined by the bulk gap closing 
mechanism [22]. The TSF resembles the SF with /-wave 
paring symmetry, which is consistent with the geomet- 
rical symmetry of the TOLs. The chiral /-wave TSF is 
fully gapped in bulk and has three chiral gapless edge 
states located on the boundary. More interestingly, the 
TSF can be modulated through initializing the lasers. 
Furthermore, there is one stable Majorana fermion bound 
to each vortex in the TSF, and the commensurability 
between the SF vortex lattice structures and the TOLs 
shows advantages to investigate the properties of Majo- 
rana fermions. Hence, these properties make the system 
a potential candidate to perform QTC. 

Modulation and model: Two-dimensional TOLs are 
created by three blue detuned laser beams, which have 
same wave-vector length but different polarizations and 
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are applied along three different directions: ±^e, 



^Cy and e^y, respectively. The total potential is given 
by V(r) = Vl cos(ki • r) with wave vectors ki^2 = 
A: and = k (0, 1). The pattern of the poten- 

tial is shown in Fig. 1(a), where the maxima form perfect 
TOLs. 

Now, we consider the trapped ultracold fermionic 
atoms having tripod- type level configuration (e.g., the 
lowest three Zeeman levels of ^Li atoms near the broad 
5- wave Feshbach resonance) [ill, [13, [23| shown in Fig. 
2(a). Three degenerate hyperfine ground states |1), |2) 
and 1 3) are coupled to an excited state |4) through spa- 
tially modulated two sets of lasers with the correspond- 
ing Rabi frequencies flv,i^ ^v,2 and Qy^s with v = a^b 
denoting two independent sets. The Rabi frequencies 
can be parameterized as l^^^i = sin6>^ cos ^^e**^^'^ ; 
^v,2 — sin^^ sin0^e*'^^'2; 1]^ 3 = cos^^e*'^^'^. = 
Xli=i Thus, the system can be described by a 

Hamiltonian, 
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^0 = - ^ l^rii^a - ^ + ^l-a- (1) 

z,q; i,j,a,(3 

Here, /i is the chemical potential. a] ^ is the cre- 
ation operator of atom on site i and in state \a) with 
a=l, 2, 3. tij is the hopping integral between site 
i and j. The Hamiltonian of the laser- atom cou- 
pling is i^/_a =5])^^ 4,4 <^i, 4 - hfty^adl^ai^^^ H.c.with 

5y the detuning to the excited state |4). Hi-a is 
diagonalized as Hi^a = ^i,v,a.(^i^y^a^i,v,a.- 4,^, a 

is the creation operator of atom on site i in eigen- 
state \Di^y^o) or Here, denote dark 

states with energy ^2,^,1/2=0 and |5i,v,Q;) denote bright 
states with energy £^,'^,3/4=^^ i v^?+^^^^- -^^^ dark 
states, |I)i,^,i)=sin(/)^e*'^^'^'^ |1) — cos (/)^e*'^^'^'2 |2) and 
|A,v,2)=cosl9^cos^^e*'^-'3'i |l) +cos6>^sin(/)^e*'^-'3'2 |2) - 
sin6>^ |3), where Sv^i^j=Sv^i — Syj. Note that \Di^y^i/2) 
have no coupling with the initial excited state |4), and 
their energies are not adjusted by the laser fields. There- 
fore, the dark states are stable under atomic spontaneous 
emission. With the adiabatic approximation [ill [23| . we 
can neglect the coupling between the dark states and 
bright states and reduce the Hamiltonian into the sub- 
space spanned by the dark states. In general, |l^i,a,i/2) 
are different from |l^i,5,i/2)- However, they can be equiv- 
alent through initializing the parameters of the lasers, 
and the equivalence attributes to the periodicity of TOLs. 
We modulate the two sets of lasers configuration il- 
lustrated in Fig. 2 (b) and (c) with the parameters 

Oa ^ (t^a = f , Sa^l ka^lX, Sa,2 -^a,l^, 

Sa,3 ka,3Z and Ob = h^lX^cp, (j)b = /C6,2^ + f , ^6,1 = 0, 

Sh,2 = 0, Sb,3 = kb^sz, where (p is an arbitrary phase 
and kv^2 and ky^s are the wave vectors of the lasers 
along the x, z axes, respectively. The wave vectors of 
the lasers are initialized to fulfil the relations: /c^,i=47r, 
/c^,2=47r/\/3, so that the commensuration with the TOLs 
is guaranteed simultaneously. Now, \Di^y^i) = ^{\l) — 

|2)), |A,^,2) = ^(|l) + |2))cosa-sina|3) for both v = a 
and b. That means the Hamiltonian ([T]) can be pro- 
jected into the subspace spanned by \Di^^) and 
with \Di^^) = \Di^y^i) and \Di^i) = |I^i,v,2) if the atoms 
are initially pumped to these dark states and they remain 
in the dark states. Then, in the dark states subspace, the 
projected Hamiltonian is 

-^eff ^ ^ l^^icT ^ ^ t'i jC^^^U(j (jf'i jCj (jf. (2) 

Here Tiia = c|^Ci,cr with a =ti denoting two pseudo- 
spins, cl ^ is the creation operator of atom on site i 
and in state |l^i,cr)- Ucr^cr'i,j=^^ ^""^ (Aa.^,^'+Ab,^,^/) dr 
Peierls phase factor from the gauge vector potential A, 
and Ay^c7C7'='i{Dv^cj\^\Dy^^>). 




FIG. 1: (Color online) (a). Uniform triangle lattices are 
formed from the maxima of the potential given by V(y) — 
Vl Y^^i cos(k^ • r). The defined lattice vectors and are 
shown. The hexangular zone encircles by the red-dashed line 
is the unit cell of triangular lattice, (b) The Brillouin zone 
(BZ) of triangular lattices. 




FIG. 2: (color online) Illustration of the light-atom interac- 
tion for generation of effective non-Abelian gauge fields and 
effective Zeeman fields, (a) The configuration of the hyper- 
fine levels of ultracold atom and three sets of laser beams 
characterized by the Rabi frequencies ^a,i, ^b,i and Qc,j with 
i — 1,2,3 and j — 1,2. The atom and the laser fields have 
interaction through the Raman-type coupling with a large 
single-photon detuning 6a/b/c- The laser beams configuration 
for Qa, and Qc are shown in (b), (c) and (d). (e) The 
relative energy levels modulated by the atom-laser couplings. 



We evaluate that the effective k^ RSOC can be simu- 
lated by the aforementioned laser fields. For convenience, 
we define the (next) nearest neighbor lattice vectors: 
{di) shown in Fig. 1 (a) with z=1...6, and set the lattice 
constant 1. For the two sets of lasers configuration illus- 
trated in Fig. 2 (b) and (c), only /7s,=/7t^=e^^^^°^^^- 

and /7d2=^d5~^*^^^°^^^^ ^^^^ nontrivial phase 

factors and other Us^/di are trivial and equal 1. Since 
the TOLs have the rotation symmetry of point group 
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C2,v , rotating the laser beams or the lattice systems with 
±^ gives another two groups of nontrivial phase fac- 
tors: /7s3=/74=e^'^^^°^^^^ /7d,=/7^^=e-^'^^^°^^^^ and i^t(k) = 



S6 

-i47r cos ^(Tx 



t 



3 if ay 



respec- 



tively. In the following section, we will show that the 
joint effects of all nontrivial phase factors Ua^a'ij are 
equivalent to the RSOC. 

In order to generate the effective ZF to split two 
pseudo-spin states, two additional laser beams that cou- 
ple the states |1) and |2) to the excited state |4) with 
a large detuning Sc are applied [24]. The laser- atom in- 

teraction is i^/_^=X](5ca{ 4ai,4 - ^^c,i<^i,4<^i,a + i^-c 

i i,a=l 

The corresponding Rabi frequencies are parameterized 



as 17c, 1 
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with 



^c,i|^ + |^c,2|^and hQc <^ The lasers config- 
uration is illustrated in Fig. 2 (d). H[_^ induces a pertur- 
bation Hamiltonian involving ground states |1) and |2), 



-hnpe-'^'''''al^ai^2 + i^.c.with Qp 



under 



the conditions of <C \sb^A\ < \£a,A\^ which guarantee 
that Hp cannot pump the atoms outside of the dark-state 
subspace. Furthermore, Hp results in a splitting between 
|A,t) and I A,;) as Hs hnp{cl^^Ci^^ - cos^ (pcl^^Ci^^). 
Note that we neglect a constant ac- Stark shift term 
—hftp{al^ai^i + aj2<^i,2), whose effect can be canceled 
with a frequency offset of the laser beams 17^,3 applied to 
the level |3). Now, the new Hamiltonian: Hq = Heff-\- 
H.. 



(3) 



Here, we renormalize the Hg and the chemical potential 
terms, /i' = /i+ hQp{l — cos^(/p)/2 and hz = hQp{l + 
cos^ (p)/2. 

Topological SF and Majorana Fermion: The SF states 
can be induced by atomic interaction from the s-wave 
scattering. The interaction term is described by the 

Hamiltonian H^ Y.i Y.a<p Va^a\^^a\^^ai^f3ai^a, where 
a and /3 label three ground states of atoms and Vaj3 are 
proportional to 5- wave scattering lengths between a, {3 
channel. In s-wave SF state. Hint can be decoupled on 
the mean-field level. Hmf = Ya<(3 
with Ac/3 = Vai3 {cii,(3Cii,a) the SF order parameter. Un- 
der the condition of Aq,/^ <C ^a/hi it is safe to consider 
the SF in the dark-state subspace, because Hmf can- 
not pump the atoms outside of the dark-state subspace. 
Then, we project Hmf to the dark-state subspace and 
have H'^j = AocJ^cJ^ + i^.c. with Aq the linear com- 
binations of Aa(3. Considering the hopping up to next 
neighbor order, the total Hamiltonian Hf = Hq 
After the Fourier transformation, in momentum Nambu 
bases: [ck^, CH,c^k^, c^j^^, ]-^, Ht can be expressed as: 



hzCTz + gk • Cr 

iAoay 



-(Sk- 



-iAoay 



Here 



ek= 



-2ti(cos A:^+2 cos ^ cos 



(4) 



2^3 (cos V^ky^2 cos ^ cos ^) 
cr={(Ja,,(jy) ak=2t2(sin/c^ - 
6k = 2t4(— sin V^/c2y+2 sin 



2 sin ^ cos ^ 



cos 



gk=(ak,^k), 

) and 
) , in which 



y I ^ 2 2 

ti=tcos(47r COS (/:)), t2=t sin (47r cos (/:)), ^3=^' cos(47r cos (^), 
t4=t' sin(47r cos with the (next) nearest neighbor 
integral t{t^). 

Now, we give the estimations about the parameters re- 
lated to the aforementioned simulations to ensure the 
experimental feasibility and rationality. For trapped 
atoms: ^Li, the wave length of laser beams utilized to 



IS aL 
kHZ with kL 



1/im, and the lattice constant 
- nx27rx30 



2m 
) 46 



produce the TOL is Xl 

. The recoil energy is Er 

1^, and t = ^(^^ 
when Vo » Er with Vp ^ W l the depth of the TOL. 

kinetic energy of atom and the renormalized factor. The 
typical atomic velocity is about several centimeters per 
second, and Ekin has the same order of ^ depends 
on geometry of the lattice. According to the ID lattice 
results [26] and 7^ < 1, we estimate that Vq ~ ?>Er is 
enough to get t' /t = Actually, the TSF is robust 

even when t' /t ^ 10~^. Here, without lost of generality, 
we set t' /t = Then t ~ 0.1 with the aforemen- 

tioned formula. From the harmonic-potential approxi- 
mation, the energies of the atoms tightly confined at a 
single lattice site are quantized to levels seperated by 

2iTfyjjQ=2Er\r^^M- On the other hand, the maximal 



band width for the TOL is Wm ^ lOt Er. It is safe to 
describe the system with single-band approximation be- 
cause of 27rhuJo > Wm- The Rabi frequencies ^a/b/c are 
10^ Er/h^ and the Qp can be tuned from to Er/h which 
is enough for Q„ <C 1^6,41 /h < \Sa,4 \ /h. Then, adiabatic 
approximation [2^ is reasonable. The typical s-wave pair- 
ing potential A^/s in experiments is about Q.l£'r-//i[27|, 
which is much smaller than Cla/b- Hence, our proposal 
is experimentally feasible when the parameters lie in the 
estimated region. 

For convenience to discuss the properties 
of SF, we rewrite ^t(k) in the new bases 
[V^k,+ , V^lk,+' -V^-k,-' -V^k with V^k,± = 



1 

V2 



Ckt ± c. 



-ki 



, Hj. has the following form: 



iey^cTy H- (k) 



(5) 



with i7±(k)=± [{hz T Aq)(7z + akcr^ ± hkCJy]- The spec- 
trums of Hamiltonian dSJ are ±£^±(/c), and 
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E±{k) = ^Jel + Al + el±2^hlAl + elel. (6) 



Here, 0^ = y^^ + ISkl^- The topological transition 
point is determined by bulk energy gap closing condi- 
tion: E-{k) = (Fig. 3(d)). The spectrums are fully 
gapped off this p oint, and the SF is topologically nontriv- 
ial when hz > y^Aq + eklk=(o,o) (Fig. 3(c)) and normal 
when hz < V'Ag + 4lk=(o,o) (Fig. 3(e)). Around the F 
point in BZ, 

n±{k) = t[(Ao T hz)cTz + Xsoikta^ + fc^a^)], (7) 

where cr± = ^ {ax ± (Jy)^ ^± = i i^x i ky). The term 
Xso{k^o-± + ^+cr^) is the RSOC with ampli- 
tude: Aso=y. It is explicit that H±(k) has the 
well-defined /-wave chirality. That means the 
topological Chern number [28] can be calculated 
Cn+ = ^ [sign(/i^ - Ao)/i,>Ao - sign(/i^ - Ao)h,<Ao]=^ 
while C'H_=0. From the square lattice results that Aq 
has maximum when the filling ^ 1 atom per site [29], 
we set Aq = 0.5t and hQp = 3Ao for the chemical 
potential around jui and filling about 0.6 atom per site. 
When the chemical potential locates at the region of /i4 
(Fig. 3 (b)), the low-energy behaviors are dominated 
by i^;(k) with k - (0,^^). Then we cannot define 
the specific chirality, and we ascribe these cases to NSF. 
According to the aforementioned analysis, we draw the 
phase diagram in Fig. 3(e). We find that the TSF 
strongly depends on the initial parameter (p of the laser 
beams and the fillings. That means we can control 
the topological properties of the SF by modulating the 
parameters of the lasers. That provides convenience to 
investigate the TSF. 

In lattice case, the ground state Chern 
number of Ht can be calculated with 
Cn=^lBzd'k2lm{^f-^), where u„(k) is 
the ground state wave-function for the nth occupied 
band (n=l,2). The straightforward calculation gives 
Ci=0 and C2=3 when hz > y^Aq + s'i\h={o,o), and 

Ci=0 and C2=0 when hz < yAfTiJ|k=(o,o)- The 
non-zero Chern number means the same numbers of 
gapless edge states from the bulk-edge correspondence. 
From the energy spectrum of Ht(ky^x) shown in 
Fig. 4, we find that three gapless chiral edge states 
transport on one edge of TSF. (see Fig. 4 (a) and 
(d)). The effective Hamiltonian describing the chi- 
ral edge states is Uedge=Y.ky>Q'^^K^kyi^)^ky{x), 
where is the effective velocity at F, and 

The Majorana condition requires ^ -ky{x) = (x), 
i.e., U-ky^cri^) — ^ky ^(^)- check that it's indeed the 
case in our model. That means three edge states are 
chiral Majorana edge states. 




FIG. 3: (color online) (a) The band structures oi Ek = £fcitGfc 
with cosip = —0.101 and /i = —2.59. (b) The band structures 
along high symmetric lines (Fig. 2 (b)). The green/dashed and 
bule/dashed lines are Sk ± \gk\ and the red lines correspond to 
(a). Four different fillings with chemical potential /ii/2/3/4 = 
-2.59t, -3.436t, -3.89t, 0.45t are shown, (c) (d) (e) are 
the SF quasi-particle spectrums corresponding to /xi, /i2 and 
/is. (f) The phase diagram as change of cosip and /i. Two 
different phases, NSF and TSF are identified. Aq = 0.5t and 
hQp — 3Ao. 




FIG. 4: (color online) The spectrums of Hamiltonian (|4]) with 
edges at x direction ix G (1,31). (a) (b) and (c) correspond 
to /ii /i2 and /is cases in Fig. 3 (b). The green/dashed lines 
indicate the contributions from the first BZ. (d) The edge 
states (the red states in (a)) transport along the edges. 



In general, the topological defects bound Majorana 
zero modes in TSC/TSF[l|, [30|. Here, we consider the 
topological excitations of vortex structures in our sys- 
tem. The quasiparticle excitations usually are described 
by the Bogoliubov-de Gennes (BdG) equation. From Fig. 
4(a), we can find that the wave function of low energy 
excitations can be constructed from the contributions of 
quasiparticle around F. For simplicity, only the non- 
trivial part l-L-\-{k) of Eq.© is taken into account. Define 
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z = X -\- iy, then d^/z* = dx ^ '^dy The BdG equation 
for the quasiparticle has the form: H+(2:,2:*)^o = £^^o 
with ^0 = ['^O'^o]^ corresponding quasiparticle 



creation operator: = y dzdz"^ '^0'^+ + '^0'^+ 
uniform TSF states, we can assume a trivial wave func- 
tion: uo = e^''^^z~ie~i'^ ^^fi = as a 
test wave function from the BdG equation. The SF order 
parameter with a vortex structure can be approximately 
expressed as Ao(r)=0 for r < Vc and Ao(r)=Aoe*^ for 
r > Vc with vorticity 1. We imagine that the vortex 
is created adiabatically by changing the wave function 
slowly enough so that it always remains an eigenstate. 
The wave function of the vortex state can be obtained 
from a singular gauge transformation: ^0 ^ ^oe*^^ with 
g = ±1 identifying the quasi-particle and quasi- anti- 
particle. In reverse, the vortex can be gauged away by 
the inverse singular gauge transformation: k k — V6/2 
and Aoe*^ Aq [31|. After the inverse transformation, 
the state is one of eigenstates, namely, the vortex ex- 
cited state. In analogy to the Laughlin's argument [32] 
about the vortex excitation in quantum Hall state, we 
get the wave function describing vortex zero mode as 



In the 



3 _ 

" 2 e 



and v'q = {u^qY. The 
unique one zero mode for /-wave case is proven by the 
numerical calculation [8]. 

The stability of Major ana zero mode is measured by 
the minigap Eg ~ A^/Ef with Ef the Fermi energy. In 



our case, the Ef is measured by t not Ekin- Hence, The 
ratio Eg /Aq can be large enough to protect Majorana 
zero mode. Take half filling as an example, we assume 
the optimized Aq ~ t and roughly estimate Ef ~ 3t. 
Then Eg /An ~ 1/3. Comparing with the hybridized 
systems p-Q, the energy scale of Ef has the order of 
electrons' kinetic energy Ekin and Aq from the proximity 
effect is much smaller compared to Ekin- So, the minigap 
in hybridized systems may be relative small compare to 
superconductive gap Aq. 

Conclusion: In summary, we have proposed a scheme 
to produce RSOC and ZF through the laser-atom in- 
teraction in TOLs, and a novelly chiral /-wave TSF is re- 
alized thanks to the 5- wave Feshbach resonance. We find 
that there exists three Majorana edge states locating on 
the boundary of the system and one Majorana fermion 
bounding to each vortex in the TSF state. The TSF 
can be controlled by modulating the parameters of the 
laser. The controllability provides convenience to inves- 
tigate the properties of the TSF. Our proposal enlarges 
TSF family and presents some advantages to study the 
Majorana fermions. 
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